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Abstract. We consider ergodic random Schrodinger operators on the metric graph Z d with 
random potentials and random boundary conditions taking values in a finite set. We show 
that normalized finite volume eigenvalue counting functions converge to a limit uniformly in 
the energy variable. This limit, the integrated density of states, can be expressed by a closed 
Shubin-Pastur type trace formula. It supports the spectrum and its points of discontinuity 
are characterized by existence of compactly supported eigenfunctions. Among other examples 
we discuss random magnetic fields and percolation models. 



1. Introduction 

This paper is devoted to the spectral analysis of certain random ergodic operators defined 
on metric graphs with a Z rf -structure. More precisely, we consider a Schrodinger operator 
with an ergodic random potential which takes values in a finite set of potentials, with ergodic 
random boundary conditions from a finite set, and with ergodic random magnetic fields from 
a finite set. Such a model may be called random quantum graph. A restriction of such an 
operator to a finite cube with selfadjoint boundary conditions has purely discrete spectrum, 
hence it has a well defined eigenvalue counting function. As the size of the cube tends 
to infinity the normalized counting functions converge to the so-called integrated density of 
states (IDS) or spectral distribution function. The corresponding result is well established for 
ergodic Schrodinger operators in the continuum (i.e. on L 2 (]R 2 )) and on the lattice (i.e. on 
£ 2 (1 d )), see for instance |Pas80l IKM821 IKir89l ICL901 IFF92] . A modification of the proof 
to ergodic, random Schrodinger on metric graphs was presented in [HV, § 6] and yields the 
existence of the IDS. 

However, the abovementioned results concern pointwise convergence of the distribution 
functions, or rather weak convergence of the corresponding measures. In the present paper 
we show that for the class of random operators at hand the convergence holds uniformly in 
the energy variable. 

Two features of the operators considered are helpful to derive such a strong statement 
about the convergence: Firstly, since the potential and other data take only finitely many 
values frequencies of finite configuration patterns exist in a nice way. Secondly, since the 
underlying geometric structure exhibits in many ways one-dimensional properties the relevant 
perturbation operators turn out to have a uniformly bounded spectral shift function. The 
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latter feature can be understood as a finite-rank property in a generalized sense, see Section 
[3] for a more precise statement. 

These two properties of the considered operators enable us to apply an abstract ergodic 
theorem for Banach space valued random variables which was proven recently in [LMV] . 
There the ergodic theorem was already used to establish uniform convergence for a certain 
class of random Schrodinger operators on £ 2 (Z d ). The situation for Schrodinger operators on 
metric graphs is somewhat more complicated since the operators are unbounded and the IDS 
is not a probability distribution any more. Still, quantum graph operators are technically 
much easier to handle than Schrodinger operators in truly multi-dimensional space, i. e. on 
L d (R 2 ),d > 2. 

The paper is organized as follows. In Section[2]we define precisely the operators we want to 
study and state the main results. In Section [3] we discuss finite rank perturbations of the kind 
as they appear in our proofs. In Section [4] we show that the ergodic theorem of [LMVj can be 
applied to the sequence of spectral shift functions obtained from an exhaustion of the metric 
graph by cubes. This is then used in Section [5] to prove our main results. In Section [6] we 
show how magnetic fields on metric graphs manifest themselves in the boundary conditions. 
Finally, in Section [7] we discuss site- and edge-percolation as an example for our results. 

2. Model and results 

We define a metric graph Gd over % d in the following way. Let ej, j = 1, . . . , d be the 
standard basis of the real d-dimensional space Mr. The vertex set Vd of the graph consists of 
the points Z d C M d which have integer coefficients. Now, to each vertex x and j G {1, . . . , d} 
we associate the edge e = [x,x + ej] with starting point s(e) = x and endpoint r(e) = x + ej 
given by 

[x, x + ej] := {x + tej : t e]0, 1[}. 

Thus, each edge can be canonically identified with the interval ]0,1[. This procedure may 
seem to induce an orientation on our graph. However, it turns out that all relevant quantities 
are independent of the choice of orientation. 

The set of all edges is denoted by Ed- The union Vd U [j eeEd e is a closed subset of M , 
hence we can consider it as a topological space with the metric inherited from M. d . We denote 
this metric subspace by Gd- We will need to consider finite subgraphs G of Gd as well. By 
a subgraph we mean a subset of the edges of Gd together with all adjacent vertices. All 
functions we consider live on the topological space Gd or subgraphs of it. This shows, in 
particular, that the choice of orientation at each edge is not relevant. 

The operators we are interested in will be defined on the Hilbert space 

L\E d ) := © L 2 (e) 

and their domains of definition will be subspaces of 

W 2 < 2 (E d ) := W 2 ' 2 (e), 

e&E d 

where W 2,2 {e) is the usual Sobolev space of L 2 (e) functions whose (weak) derivatives up to 
order two are in L 2 (e) as well. The restriction of / £ W 2,2 (Ed) to an edge e is denoted by f e . 
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For an edge e = [s(e), r(e)] = [x, x + ej] and g in W 2,2 (e) the boundary values of g 
g{s(e)) := ^9(x + tej), g{r{e)) := limg(s(e) + tej) 

and the boundary values of g' 

5 (s(e)J := hm - and g (r(e)) := lim — 

e\o e e\o — e 

exist by standard Sobolev type theorems. Note that we have introduced a sign in the derivative 
at the endpoint of an edge. This makes our definition of derivative canonical, i.e. independent 
of the choice of orientation of the edge. For / G W 2,2 (Ed) and each vertex x we gather the 
boundary values of f e {x) over all edges e adjacent to x in a vector f(x). Similarly, we gather 
the boundary values of f' e {x) over all edges e adjacent to x in a vector f'(x). 

Given the boundary values of functions, we can now dicuss the concept of boundary con- 
dition. Here we use material from [KS99, I larOO to which we refer for further details and 
proofs. A single-vertex boundary condition at x G V is a choice of subspace L x of C 4rf with 
dimension 2d such that 

T]((v, v'), (w, w')) := (v, w) - {v, w') 

vanishes for all (v, v'), (w, w') G L x . An / G W 2,2 (Ed) is said to satisfy the single-vertex 
boundary condition L x at x if (f(x),f'(x)) belongs to L x . A field of single-vertex boundary 
conditions L := {L x : x £ Vj} will be called boundary condition. Given such a field, we 
obtain a selfadjoint realization of the Laplacian A on L 2 (Eii) by choosing the domain 

V(A L ) := {/ G W 2 > 2 (E d ) : Vx : (f(x), f{x)) G L x }. 

Particularly relevant boundary conditions are Dirichlet boundary conditions with subspace 
L D consisting of all those (v , v ') with v = 0, Neumann conditions with subspace L N consisting 
of all those (v, v') with v' = 0, and Kirchhoff (also known as free) boundary conditions L K 
consisting of all (v, v') with v having all components equal and v' having the sum over its 
components equal to 0. 

Remark 1. More general types of boundary conditions are conceivable. We will restrict 
ourselves to those just introduced. They may be called graph local boundary conditions, as 
each boundary condition involves values of / at one vertex only. 

Everything discussed so far including existence of limits of functions at the vertices and 
the notions of boundary condition extends in the obvious way to subgraphs. Moreover, for a 
subgraph G of Ga with edge set E, we write W 2,2 {E) := ® e& EW 2 ' 2 (e). The number of edges 
of a finite subgraph G of Gd is denoted by \E\. 

In order to define random operators we need some further data including a probability 
space (fi,P) and an action of (a subgroup of) the automorphism group of Gd on Q and maps 
L, V from Q, into the space of boundary conditions and potentials, respectively. 

For us two groups will be relevant, the full automorphism group and the group r^ rans 
of translations by Z, d . Note that ^f u \\ is generated by translations by vectors in Z, d and a 
finite set of rotations. We choose T to be ^f u \\ or r^ rans and assume that it acts ergodically 
on (O, IP) via measure preserving transformations. To simplify the notation we identify 7 G T 
with the associated measure preserving transformation. 

Let us describe the type of random operators we will consider in this paper: 
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Assumption 1. Let (f^,P) be a probability space and T G {^fulb ^transi a g rou P acting 
ergodically on P). Let B be a finite subset of -L°°(0, 1) and C a finite set of boundary 
conditions. A random potential is a map 



V : n — > B with y( 7 (w)) 7(e ) = F(w) e 

for all 7 G r and e £ E^. A random boundary condition is a map 

L:ft — with L( 7 (o;))(7(x)) = £(a/)(a;) 

for all 7 G T and u G V^. 

A family of random operators (H u ) on L 2 (Ed) can be defined with domain of definition 

V(H U ) := {/ G W 2 ' 2 (£ d ) : (f(x),f'(x)) G for all x G F d } 

acting by 

(H u f)(e) :=-f: + V[u) e f e 
for each edge e. These are selfadjoint lower bounded operators. 

We assume throughout the paper that Assumption Q] holds, and for this reason do not 
repeat it in every lemma. Note however, that in many statements we will only need part of 
the structure described in Assumption [TJ 

Remark 2. While ^( u \\ is not commutative it is a natural object to deal with. In particular, let 
us note that the Laplacian without potential and boundary conditions in all vertices identical 
to Kirchhoff conditions is invariant under Lf u jj. 

We will need to consider restrictions of our operators to finite subgraphs. These are finite 
subgraphs associated to finite subsets of Z. The cardinality of a finite subset Q of Z d is 
denoted by \Q\. Let 

d 

S:={^2t jer .0<t,<l}. 
i=i 

For a finite subset Q of Z rf , we define the associated subgraph Gq of Gd by 

G Q := (Q + S)n G d . 

The vertices and edges of Gd contained in Gq are denoted by Vq and Eq, respectively. Note 
that Vq D Q. The set VA of inner vertices of Gq is then given by those vertices of Gq all of 
whose adjacent edges (in Gd) are contained in Gq. The set of inner edges Eq of Gq is given 
by those edges whose both endpoints are inner. The vertices of Gq which are not inner are 
called boundary vertices. The set of all boundary vertices is denoted by Vq. Similarly, the 
set of edges which are not inner, is denoted by Eq. 

The restriction H$ of the random operator to Gq has domain given by 

V(HS) ■= {/ G W 2 ' 2 (E Q ) |Vx G V% : (/(*),/'(*)) G L(co)(x), 

VxeV°:(f(x)J'(x))£L D } 

This operator is again selfadjoint, lower bounded, and has purely discrete spectrum. Let us 
enumerate the eigenvalues of H$ in ascending order 

\ 1 (hS)<x 2 (hS)<x 3 (hS)<... 
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and counting multiplicities. Then, the eigenvalue counting function rQ on R defined by 

n£(A) := H{n G N | \ n (H®) < A} 

is monotone increasing and right continuous, i.e. a distribution function, which is associated 
to a pure point measure, [Q. Denote by 

N${\) := T^-rTig(A) 
\ h Q\ 

the volume-scaled version of rC(A) and note that \Eq\ = d\Q\ as the edge to vertex ratio in 
the graph (Vd, Ed) is equal to d. 

\v a I 

A sequence {Qi)i&\ of finite subsets of Z rf is called a van Hove sequence in Z d if liim_ i . 00 -i^r = 
0. For a finite subgraph H of Gd let Xff the multiplication operator by the characteristic func- 
tion of H. Denote the trace on the operators on l?{Ed) by Tr[. . .]. We can now state the 
main result of the paper. 

Theorem 3. Let Q be a finite subset of TL d . Then, the function N defined by 

(1) N{\) : = J- /" Tr [xg q X]-oc,x]{H u )] dF(co) 

I Q I J f2 

does not depend on the choice of Q, is the distribution function of a measure fi, and for any 
van Hove sequence (Qi) in 1i d 

lim ||iV* - JVHoo = 

l— >oo 

for almost every lo G O. 7n particular, for almost every lo £ Q., N$ l (\) converges as I — > oo 
pointwise to N(X) for every A G M. 

Remark 4. (a) We chose to define the integrated density of states by a Shubin-Pastur formula 
and to prove that it coincides with the almost sure limit along van Hove sequences using an 
ergodic theorem. Alternatively, one could take the point of view that the intuitively relevant 
objects are the normalised eigenvalue counting functions on finite graphs. Then the IDS 
would be defined as their limit along a van Hove sequence and the trace formula ([T]) would 
be the result of the theorem. 

(b) For random Schrodinger operators with more general potentials and Kirchhoff bound- 
ary conditions on all vertices the IDS was constructed in [HV| as the pointwise almost every- 
where limit of the sequence nS 1 ■ 

(c) In the case that Y is the full automorphism group the equality 

(2) N{\) = ^L^J^Tr [xGX]-oo,A](^)] dP(u) 

holds for any finite subgraph G with edge set E{G). To see this note that by linearity 

/ Tr [xsX]-oo,A]0Hu,)] dF{u) = V / Tr \x[o,e<]X\-oa,x](.H u j\ dF(uj) 
Jn ~[ Jn L J 

= d / Tr [x[o,ei]X]-oo ) A](- H w)] dF(uj). 
Jn 

For the second equality we used that the d terms in the sum are all equal, since the group 
r^ull acts transitively on the edges. For an arbitrary finite graph G one can do a similar 
calculation, and formula (O) follows. 
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While the definition of the IDS involves an ergodic theorem, there are other spectral features 
of H u whose almost sure independence of to uses only the ergodicity of the group action. 
Prominent examples are the spectrum cr(i? w ) and its subsets cr pp (i/^), c sc (ff a; ), cr ac (-£^), 
^disc(Huj), Oess^Hu) according to the spectral type. In fact, by applying the general framework 
of |LPV07| we immediately infer the following theorem. 

Theorem 5. There exist subsets of the real line X, T, pp , X sc , X ac , X^c, T, ess and an £1' C fl 
of full measure such that ^(H^) = X and a^H^) = X. for all these spectral types • G 
{pp, sc, ac, disc, ess} and all wgO'. 

Our goal here is to establish a relation between the almost sure spectrum of H u and its 
IDS. More precisely, in two corollaries to Theorem [3] we relate the topological support supp/U 
of \x and the set S p (fi) := {A G R : /i({A}) > 0} of atoms of fi to the spectrum of H^. Note 
that the set of discontinuities of the IDS is precisely S p (fj,). 

Corollary 6. X equals the topological support supp/i of \x. 

As usual an / G W 2,2 (E c i) is said to be compactly supported if f e = for all but finitely 
many edges e. 

Corollary 7. Denote by ^ C mp the set of energies A G R such that there exists almost surely 
a compactly supported L 2 -eigenf unction f w with H^f^ = A/ w . Then 

(3) Sp(fi) = Yj cm p. 

Remark 8. (a) We will prove equation ([3]) using the following ingredients: Theorem [3] and the 
reasoning presented in [KLS03] , Note, however, that both of the sets in ([3]) are given in terms 
of the infinite- volume operator H w rather than by its restrictions to finite subgraphs. Indeed, 
there exists a proof of ([3|) which makes use neither of finite-graph operators, respectively 
their eigenvalue counting functions, nor of an ergodic theorem: Once one has established the 
invariance property of the IDS, as done in e.g. [LP V07] . equation ([3j) follows from the mere 
existence of a van Hove sequence. This strategy is implemented in the proof of Theorem 
2.3 (i) of [Ves05b| for operators on combinatorial graphs and can also be adapted to our 
situation. Note however, that there are certain situations where it is rather natural to relate 
the (size o f the) jumps of the IDS to those of its finite-graph analogues, see for instance 
|DLM+03llElellL\7] . 

(b) There are many examples where the IDS has discontinuities. Indeed, the free Laplacian, 
i.e. the Schrodinger operator with identically vanishing potential, with Dirichlet, Neumann, or 
Kirchhoff boundary conditions exhibits compactly supported eigenfunctions. The relevance 
of Theorem [3] is that one has uniform convergence all the same. 

(c) If the randomness entering the potential of the operator is sufficiently strong it is 
natural to expect a smoothing effect on the IDS. In fact, in [HV| for a class of alloy-type 
random potentials the Lipschitz-continuity of the IDS was established. In [GV] we show for 
a different class of random potentials how one can estimate the modulus of continuity of the 
IDS. 

3. Finite rank perturbations 

We will have to deal with changes in the spectral counting function due to changes of 
boundary conditions and changes of potentials. In particular, we will have to cut graphs 
into smaller pieces by suitable boundary conditions in order to apply the ergodic theorem of 
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[LMVJ . These changes are finite rank perturbations in an appropriate sense and hence change 
the spectral counting functions only in a very quantifiable way, which is made precise in this 
section. 

Let Hi and H 2 be selfadjoint lower bounded operators with discrete spectra. Thus the 
corresponding counting functions uh 1 and njj 2 are well-defined. Then, the spectral shift 
function is the difference of two counting functions, i.e. 

We will make use of the min-max principle in the form 

n H (X) = max{dimX | X < V(H),H\ X < A} 

where X runs over all linear subspaces (denoted by <) of the domain of H on which H is 
bounded above by A in the sense that (/, Hf) < A(/, /) for all / £ X. The following lemma 
will be useful for dealing with several types of perturbations: 

Lemma 9. Let Hi, H2 be selfadjoint lower bounded operators with discrete spectra on a 
separable Hilbert space. Assume that T>q < T>(Hi) PiT>(H 2 ) such that Dq has finite index k in 
V{H 2 ), i.e. dimV(H 2 )/V = k. If H 2 \ Vo > H\\v then 

nu 2 W < "Hi (A) + k, i.e. (A) < k. 

Proof. By the min-max principle for the counting functions we have: 

n H2 (X) = max{dimX | X < V(H 2 ),H 2 \ X < A} 

< max{dimX | X < T> , H 2 \ x < A} + k 

< max{dimX | X < V ,Hi\ x < A} + k 

< max{dimX | X < V{Hx),Hi\ x < A} + k 
= n Hl (X) + k 

and the proof is finished. □ 

This easily leads to the familiar boundedness of the spectral shift function for finite rank 
perturbations: 

Corollary 10. Let H be selfadjoint with discrete spectrum and K selfadjoint and bounded. Let 
Y + (resp. Y-) be the strictly positive (resp. negative) spectral subspace of K. Ifk + := dimy + < 
00 then n(H, X)-k+ < n(H + K,X). Ifk--.= dim F_ < 00 then n{H + K,X) < n(H, A) + /c_. 
In particular, the spectral shift function for finite rank perturbations is bounded by the rank 
of the perturbation. 

Proof. Applying the previous lemma to H 2 = H + K, Hi = H (such that D(H 2 ) = T>{H{) = 
T>(Hi)), choosing T>q = T>{H) n Yjt, yields one half of the claimed estimate. The other half 
follows if we replace K by — K and H by H + K. □ 

The perturbations in our applications will be of finite rank only in the following general 
sense: 

Corollary 11. Let Hq be densely defined, closed, symmetric and lower bounded with defi- 
ciency index k. Let Hi, H 2 be two selfadjoint extensions of Ho with discrete spectrum. Then, 

(4) ieffi^ 2 |<fe- 

Proof. This is a direct consequence of the previous lemma with Dq = D{H$). □ 
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Remark 12. Alternatively this corollary could be proven using the invariance principle for 
the spectral shift function and the second resolvent identity. In fact, this would not need the 
assumption of discrete spectra. 

Let us now turn to the operators of interest for our applications. These are the restrictions 
of —A to a finite subgraph of Gd, with selfadjoint boundary conditions and possibly with an 
additional L°°-potential. 

Lemma 13. Let G be a subgraph of Gd- Let A and A' be two selfadjoint realizations of the 
Laplacian on G. Set T> = T>(-A) n V(-A') and k := dim£>(-A)/£> - Then k < 2\E\, and 
|f-A,-A'| < k 

Proof. Set m := \E\. Obviously W^ 2 {E) := {/ G W 2 ' 2 (E) : (f(x),f(x)) =0ViG V(G)} 
is a subspace contained in all possible domains, and all of them are contained in W 2,2 (E). 
Because of dimW 2 ' 2 (E) /W ' (E) = Am and the rank 2m difference for the domains, T>$ has 
index at most 2m, and we can apply Lemma [HJ □ 

We are particularly interested in Dirichlet resp. Neumann conditions on the edges of a 
subgraph G. The corresponding operators — A^, — A^ decompose completely as direct sums 
over the edges, and the counting functions are given by ra_ A g(A) = \E\n£)(X) and n_ A G(X) = 
\E\nj^(X). Here, 

A , (Lv^/ttJ if A > 0, 
n D{X) = < n N {X) = n D (X) + X[o,oo[(A) 

I (J otherwise, 

are the counting functions for the Dirichlet and Neumann Laplacians on ]0, 1[. [x\ is defined 
for x £ Mby [x\ £ Z and x — 1 < [^J < x. Together with the previous lemma, this implies 
the following proposition. 

Proposition 14. Let G be a finite subgraph of Gd and A any selfadjoint realization of the 
Laplacian on G. Then, \ri-^(X) — \E\ud{X)\ < 2\E\. 

In fact, looking at the quadratic forms one sees that the Dirichlet form dominates the forms 
of all other selfadjoint extensions; therefore, \E\nn is a lower bound to the counting functions 
of all other selfadjoint extensions. 

Of course, everything said so far about selfadjoint extensions remains true when we add 
bounded perturbations. The spectral shift is determined in the following lemma. 

Lemma 15. Let —A be a selfadjoint realization of the Laplacian on a subgraph G of Gd and 
let W\,W2 be bounded selfadjoint operators on L 2 (E). Set Hj = —A + Wj,j = 1,2. Then 
there is a constant C depending only on the norms ofW\,W2 such that 

\Z Hl ,H 2 \(X)<C\E\ 

Proof. Set m = \E\. According to Proposition 1141 

m( — 2 + n£)(X)) < n_A(A) <m(2 + no(A)) and therefore 
m(-2 + n D (X-\\W\\)) < n. A+w (X) < m(2 + n D (X + \\W\\)) 

for any bounded selfadjoint W because — ||W|| < W < \\W\\ implies n#(A — ||W||) < 
riH+w(X) < uh{X + \\W\\) for any selfadjoint H. Consequently, the explicit form of no 
given above leads to 

m( - 3 - vWliA) < «-A+w(A) - v^/tt <m(2 + ^W\\/^) 
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and subtracting the estimates for Hi and H2 we finally get 

\t Hl ,H 2 \ (A) < m(5 + VWJA + y/W^M 

Using the fact that the Dirichlet counting function is a lower bound would allow us to replace 
5 by 3. □ 

Remark 16. The considerations of this section hold for general finite metric graphs with all 
edges of equal length. 

4. An abstract ergodic theorem 

In this section, we discuss the abstract ergodic theorem of |LMV] and show how it can be 
applied to give convergence of a suitably defined spectral shift function of a family of random 
operators. See [Len02| for earlier one-dimensional results of the same type. This is the crucial 
step in the proof of our main result. 

Let A be a finite set and choose T G {rtrans» ^fullJ -- Then, we can proceed as follows. The 
set of all finite subsets of 7L d is denoted by T . A map V: 7L d — > A is called an A- colouring 
of Z d . A map P: Q{P) — ► A with Q(P) G T is called an A-pattern. For MeNwe denote 
by Cm the cube at the origin with sidelength M — 1, i.e. 

C M := {x £ Z d : < xj < M - l,j = 1, . . . ,d}. 

The set of patterns P with Q(P) = C M is denoted by Vq(M). The set of all C M , M G N, 
is denoted by C, and a pattern with Q{P) G C is called a cube pattern. For a pattern P and 
Q G T with Q C Q{P) we define the restriction P n Q of P to Q in the obvious way by 
P n Q : Q — ► A, x P(x). For a pattern P and 7 G T we define the shifted map 7P by 
7-P : r yQ(P) — * A, 7(2/) 1— > P(y)- On the set of all patterns we define an equivalence relation 
by P ~ P' if and only if there exists a 7 G T with 7P = P'. For a cube pattern P and an 
arbitrary pattern P' we define the number of occurrences of the pattern P in P' by 

f P P' := G Q(P') : P' n (x + Q(P)) ~ P}. 

Definition 17. A map 6: J 7 — ► [0, 00) is called a boundary term if = 6(i + Q) for all 
i G Z d and Qef, lim^oo |Qj = for any van Hove sequence (Qj), and there exists 

P> > with 6(Q) < D|Q| for all Q £ J 7 . 

Definition 18. Let (X, \\ ■ ||) be a Banach space and F: T — + X be given. Let V G 

{ r trans> r full} be S iven - 

(a) The function F is said to be almost-additive if there exists a boundary term b such that 

m 



k=i 



F(U% =1 Q k )-Y,F(Qk) 

fe=i 

for all m G N and all pairwise disjoint sets Qk G k = 1, . . . , m. 

(b) Let V : 7L d — > be a colouring. The function F is said to be T -V -invariant if 

F(Q) = P( 7 Q) 

whenever 7 G T and Q £ J 7 obey 7(V n Q) = V n (7Q). In this case there exists a function P 
on the cubes C with values in X such that 



P( 7 Q) = p( 7 - 1 (Vn( 7 Q) 
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for cubes Q £ C and 7 £ T. 

(c) The function F is said to be bounded if there exists a finite constant C > such that 

<C|Q| 

for all Q £ J 7 . 

Theorem 19. Let A be a finite set, V: Z d — ► A an A-colouring and (X, \\ • ||) a Banach- 
space. Let T £ {I *f u n,I trans} ^ e given. Let (Qj)jen be a van Hove sequence such that for 
every pattern P the frequency vp = linij^oo |Qj| _1 t)p(V Pi Qj) exists. Let F : T — > X be a 
T-V -invariant, almost- additive bounded function. Then the limits 

j-» m «-» pe ^ M) \C M \ 

exist and are equal. 

Remark 20. (a) The theorem is proven in [LMV| for V = r^ rans , However, the proof carries 
over verbatim to V = Tf u jj. 

(b) We have explicit bounds on speed of convergence in terms of speed of convergence of 
the frequencies. For details we refer to |LMV| . 

Now we turn to the case where the colouring is given by random potentials and boundary 
conditions. Then each induces a colouring V(u>) with values in A := (®j =1 B) x C given 

by 

V(u)(x) := {{V{[x,x + ej\)jM. 

Lemma 21. Let (Qj) be an arbitrary van Hove sequence. Then, for almost every u £ Q, the 
frequency up = lirrij^oo \Qj\~ l $ p (y(u>) fl Qj) exists and is independent of uj for every cube 
pattern P. 

Proof. For a fixed pattern P the frequency exists for almost every w by a standard ergodic 
theorem. As there are only countably many P the statement follows. □ 

In our setting we also have a family of random operators (H w ) along with the Dirichlet 
Laplacian Ap. As discussed above, these operators can be restricted to the subgraphs Gq 
induced by finite sets Q of Z rf . This yields the operator Hj and with spectral counting 
functions n2 and respectively. Now, decomposes as a direct sum of operators. Thus, 
denoting the eigenvalue counting function of the negative Dirichlet Laplacian on ]0, 1[ by 
tid(A) as above, we have = \Eq\ud = d\Q\no- The associated spectral shift function is 
given as 

£?(A) ■= r$W ~d\Q\ n D (X) = d\Q\ (iV^A) - n D (A)). 

The crucial point is that £ w falls into the framework of almost additive F introduced above. 
This is shown next. 

Lemma 22. Let (T>, \\ ■ ||oo) be the Banach space of right continuous bounded functions on 
R. Then, for each uj £ Q, the function ^ : T — * T>, Q ^ £3 , is a bounded, V(oj) invariant 
almost additive function. 
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Proof. Almost additivity: Obviously, Q i— ► n"^ is almost additive with boundary term equal 
to zero, as Dirichlet boundary conditions decouple everything. As for n u , we note that 
Q = U^L 1 (5fc with pairwise disjoint Qk, k = 1, . . . , m, induces a decomposition of Gq into GQ k , 
k = 1, . . . , m. The restriction hS of to Gq differs from the direct sum ®™ =1 H UJ Q k only by 
boundary conditions on the boundary of the GQ k , k = 1,... ,m. These boundary conditions 
yield boundary terms in the corresponding eigenvalues counting functions by Lemma [I3j 
Thus, fi w is almost additive as well. Hence, £ w is almost additive as it is a difference of almost 
additive functions. 

V{u) -invariance: This is clear from the definitions. 

Boundedness: By Lemma \13\ changing all boundary conditions to Dirichlet conditions 
introduces an error term of the order 2\E\. On the other hand, the operator with potential 
and only Dirichlet conditions can easily be compared to the Dirichlet operator without any 
potential by Lemma [T5l 

□ 

Remark 23. The need to use a spectral shift function, i.e. the difference between n w and no, 
in the above proof comes exclusively from the last step in the proof. This use of a spectral 
shift function to define the IDS bears some similarity to how the integrated density of surface 
states is defined, see for instance [EKSS881 IEKSS9U1 ICTra99al ICha99bl iKSOl] . 

The key result is now the following proposition. 

Proposition 24. There is a bounded right continuous function S: R — > R such that for a 
given van Hove sequence (Qi) for almost every the uniform convergence 



lim 

l— >oo 



r 



E c 







holds. 



Proof. Almost sure existence of the limit is a direct consequence of the previous two lemmas 
and the first theorem of this section. In fact, this theorem gives an explicit formula for the 
limit in terms of frequencies vp and F. This formula shows that the limit does not depend 
on uj. □ 



5. Proofs of the main results 

In this section we gather the material of the previous sections in order to prove Theorem 
[3] and its two corollaries. 



Proof of Theorem^ We only consider the case V = r^ rans . The case Y = T£ u ^ is similar 
and even simpler. 

We first show independence of N of the choice of Q: Set G x = x + S. The invariance 
assumption on the random operators, the invariance of the trace under unitary conjugation 
and the invariance of P under translations easily show that the expression 

Jn 

does not depend on x G Z d . Now, the claimed independence of Q follows easily. 
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We now show convergence of the \Eq\ 1 n® 1 , I G N. Proposition 1241 yields almost sure 
convergence of 

1 .&(.)= *rfi(.)-n D (.) 



\E Ql r yj \E Ql \ 

to a limit E with respect to the supremum norm. As the subtracted term no(') does not 
depend on I we obtain the desired convergence of ,i \ (~3(-) + hd- The limit is N := E + no- 

Finally, we have to show that N = N. To do so it suffices to show that the measures 
associated to ./V and N respectively are equal. This in turn follows once we show vague 
convergence. By a variant of the Stone/Weierstrass theorem it suffices therefore to show that 

(*) r^tr( XGQ; /(^)-/(^))^0, l^oo, 
I Qi I 

for all / of the form f[t) = (t — z)^ 1 for z with non-vanishing imaginary part. For Q = Qi we 
can split the graph into the two components Gq and Gj\Gq. Then, H w and H Gq '®H G 
differ only by boundary conditions on the set Vq. Thus, by the second resolvent identity 

D:=f(H u )-f(H° Q ®H^ GQ ) 

is an operator of rank at most 4d|Vg|. Moreover, D is obviously bounded by 2\Im(z)\~ 1 . 
This gives 



\tr( XGQi f(H w ) - f(HS'))\ = \tr( X G Ql (/(F w ) - f(H^ H^ Gq )))\ 



< 8d———\v£, 



\Im(z) 

As (Qi) is van Hove we obtain (*). □ 

Proof of Corollary^ This is a variant of the proof of Theorem 2 in [KLS03]. 

Let A be a point of discontinuity of /i. Thus, 6 := /i({A}) > 0. By the uniform 
convergence proven in Theorem [3] this means that the multiplicity of the eigenvalue A of hS 1 
is at least 5/2\Eq \ for I large. On the other hand, within the eigenspace of HS 1 to the 
eigenvalue A the space of / with (f(v),f'(v)) = for all v G Vq has codimension at most 
2c?|VqJ. As (Qi) is van Hove, we have 

2d\V^\<5/2\E Ql \ 

for I large enough and we obtain a compactly supported eigenfunction of H$ l to A all of 
whose boundary values are 0. Thus, this function can be extended by to give a compactly 
supported eigenfunction. 

<£^: Choose the van Hove sequence to consist of cubes centered around the origin with 
increasing sidelengths. For I G N, define the function \i on ^ by Xli^) = 1 if there exists 
a compactly supported eigenfunction of H u to A which is supported in Qi and vanishes 
identically on all edges e G Eq i and x/( w ) = otherwise. Then 

< S t := ( Xl(uW(u) 
Jn 
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for all I large enough. Fix such an I and set % := Xli & := $l an d R to be the diameter of Q[. 
Then, the ergodic theorem yields 

almost surely. Choose k$ with t^-t Yl X £Q k x(T x uj) > 5/2 for k > &q. For each k > k$ we 
consider a maximal system Mr of points in such that x(^ Xw ) = 1 an d points in have 
distance at least 2R. A direct combinatorial argument then shows that 

By construction different points in yield orthogonal compactly supported eigenfunctions 
of H u . This yields 

for k > ko and the statement follows from the uniform convergence. □ 

Proof of Corollary Due to ergodicity and the Pastur-Shubin trace formula the corollary 
follows by standard arguments. 

□ 

6. Operators with magnetic fields 

Our setup is general enough to include magnetic fields as well. Since this is not quite 
obvious from the definitions and since there is a special interest in these operators we devote 
a section to the question how magnetic fields enter through boundary conditions. 

To this end, let G be a metric graph and L a choice of boundary conditions as in Sections 
[2] and El The most general symmetric first order perturbation of — ^ on an edge e G E is, 
up to zeroth order terms, given by 

H(a) e := - - %a t 

for arbitrary real valued a e G C 1 (e), where e is the closure of the edge e, i.e. identified with 
the closed interval [0, 1]. If G is embedded, e.g. in W 1 such as in our case, it is natural to 

fix a so called vector potential A G ((C 1 (IR ci )) d and to define a e to be the component Aj of 
A if e is embedded with euclidean direction ej. This corresponds to a magnetic field 2-form 
B = Y^ij dxi A dxj on R d . We will work in the general setting, though. 

The selfadjoint realization of H (a) corresponding to L is then given by the domain 

V{H L (a)) = {/ G W 2 > 2 {E) : Vx G V : (f(x), f'(x) - i{af)(x)) G L x } 

as the usual partial integration argument shows; i.e. one has to specify mixed Dirichlet and 
(magnetic) Neumann boundary conditions as expected. 

Now define the phase along the edge e by (p e (t) = L a e (s) ds. Since each a e is real valued, 
multiplication by e lVe on each edge e defines a unitary operator U on L 2 (E) which maps 
W 2 ' 2 (E) into itself. A simple calculation shows that 

(U*H(a)U) e = 
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on W 2,2 (E). Therefore, Hi(a) e is unitarily equivalent to H e = — ^ with domain 

V{H L ) = {/ G W 2,2 (E) : Uf G V{H L (a))} 

= {/ G W 2 ' 2 (£) : Vx G V : ((£//)(*), - %{aUf){x)) G 

for some graph local boundary condition L since multiplication by U preserves graph locality. 
In order to determine L we calculate, using the notation from Section [21 

(Uf)Me)) = fe(s(e)) 

(Uf) e (r(e)) = e^ r ^f e (r(e)) 
{Uf)' e {s{e))-t{aUf){s{e))=f'Me)) 

(UfYMe)) ~ i(aUf)(r(e)) = e^ e )^(r( e )) 

Hence we define a unitary diagonal matrix u on C 4rf by setting a diagonal entry to e l< ^ s ( e " = 1 
if it corresponds to a boundary value at some s(e) and to e lip v( e " if it corresponds to a 
boundary value at some r(e). Then L = u*L. 

Note that Dirichlet and (magnetic) Neumann conditions for the magnetic operator H{a) 
transform under u into Dirichlet and (non-magnetic) Neumann conditions for H. This is 
to be expected because these boundary conditions decouple, and any magnetic field can be 
gauged away in strictly onedimensional systems. In contrast, (magnetic) Kirchhoff conditions 
for H{a) transform into completely different boundary conditions: At a vertex x, only the 
values at outgoing edges (those with s(e) = x) coincide; they, in turn, coincide with the 
values e l< ^ r ( e ^ / e (r(e)) at the incoming edges so that / need not be continuous at the vertex. 
Similarly, the derivative part of the Kirchhoff condition is modified to a weighted condition. 

Finally, let us note that the results of this section imply that the results of section [2] for 
Schrodinger operators with random potentials and random (or fixed) boundary conditions 
lead to the same results for magnetic Schrodinger operators with random potentials and 
random (or fixed) magnetic fields, specified through the phases ip e (r(e)) at the endpoints. 



7. Percolation on Gd 

In this section we discuss two models which can be seen as versions of site and edge perco- 
lation on Gd, respectively. A third model corresponds to percolation on the junctions between 
vertices and edges. These models are based solely on random boundary conditions. The po- 
tential of the operators is identically equal to zero. Unlike in the percolation models on combi- 
natorial graphs "deleted" edges are not removed completely from the graph but only cut off by 
Dirichlet boundary conditions. The reason is that removing edges would mean to remove infi- 
nite dimensional subspaces from our Hilbert space. This would result in a spectral distribution 
function which is not comparable to the one of the Laplacian with the concerned edge included. 
Hamilton operators on percolation subgraphs of combinatorial graphs have been considered in 
the literature in theoretical physics [dGLM59l IKE72j ICCF+86] . computational physics |KB02] 
(and references therein), and mathematical physics [Ves05a, Ves05b, KM06, MS07, Ves06]. 

Before giving details we would like to emphasize the following: The examples below include 
cases in which the graphs contain many finite components giving rise to compactly supported 
eigenfunctions. In particular, the integrated density of states has many discontinuities. In 
fact, in the subcritial phase the IDS is a step function. However, despite all these jumps our 
result on uniform convergence does hold! 



UNIFORM EXISTENCE OF THE IDS 



15 



We first discuss a site percolation model. The percolation process is denned by the following 
procedure: toss a (possibly biased) coin at each vertex and - according to the outcome - put 
either a Dirichlet or a Kirchhoff boundary condition on this vertex. Do this at every vertex 
independently of all the others. To be more precise, let p G (0, 1) and q = 1 — p be given. 
Let A := {L D ,L K } and the probability measure v := p8 L K + (1 — p)S l d on A be given. 
Define fl as the cartesian product space x x& y d A with product measure F := ® x &f d v - Let L 
be the stochastic process with coordinate maps L{oS)(x) := uj{x). These data yield a family 
of random operators — A w := acting like the free Laplacian with domain given by 

D(A W ) = {/ G W 2,2 (E) : (f(x)J'(x)) G <j(x) Vx G V d }. 

Intuitively, placing a Dirichlet boundary condition at a vertex means "removing" it from the 
metric graph. The 2d formerly adjacent edges have now "loose ends". A fundamental result of 
percolation theory tells us that for sufficiently small values of p the percolation graph consists 
entirely of finite components almost surely. For these values of p our Laplace operators 
decouple completely into sums of operators of the form — A G for finite connected subgraphs 
G of Gd- Here, A G acts like the free Laplacian and has Dirichlet boundary conditions on 
its deleted vertices (boundary vertices) and Kirchhoff boundary conditions in its vertices 
which have not been deleted by the percolation process (interior vertices). We introduce an 
equivalence relation on the set of connected subgraphs of Gd with a finite number of edges 
by setting G 1 ~ G 2 iff there exists a 7 G T such that 7G 1 = G 2 . For such an equivalence 
class Q we define rfi as the eigenvalue counting function of — A G for some G G G, and set 
= p^TT- Defining the density vg of an equivalence class of finite subgraphs of Gd within 
the configuration uj in the obvious way, we obtain as integrated density of states for the family 

g 

where the sum runs over all equivalence classes Q of finite connected subgraphs of Gd- Thus, 
the integrated density of states is a pure point measure in this case with many jumps. More 
interestingly, all these jumps remain present (even if their height is diminished) when we start 
increasing p. This yields models in which the operators are not given as a direct sum of finite 
graph operators but still have lots of jumps in their integrated density of states. Related 
phenomena for combinatorial Laplacians have been studied e.g. in [CCF + 86l IVes05b~l IVes06] . 

We now discuss an edge percolation model. The basic idea is to decide for each edge 
independently whether Dirichlet boundary conditions are put on both ends or not. All other 
boundary conditions are Kirchhoff type. The problem when defining this edge percolation 
model is that our stochastic processes are indexed by vertices rather than edges. We thus 
have to relate edges to vertices. This is done by going to each vertex and then tossing a 
(biased) coin for each j = 1, . . . , d to decide how to deal with the edge [x, x + ej\. 

Here are the details. Let po G (0, 1) and p% = 1 — po be given. Let A consist of all 
maps S from {1, . . . , d} to {0, 1}. Put a probability measure v on A by associating the value 
Y[j=iPs(j) to the element S. Now, Q is the cartesian product space x xe v d A with product 
measure P := ® x ^v d v - To each u G f2 we associate the operator — A w = H w which acts 
like the free Laplacian and has boundary conditions as follows: The edge e = [x, x + ej] has 
Dirichlet boundary conditions on both ends if the random variable associated to the vertex 
x has the j-th component equal to 1. Otherwise the boundary condition is chosen to be 
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Kirchhoff. Here, again the operator decouples completely into operators on finite clusters for 
small enough values of pq. 

Similary one can consider a percolation process indexed by pairs (x, e) of adjacent ver- 
tices and edges. As in the last model consider a colouring A consisting of all maps S from 
{1, —1, 2, —2 . . . , d, —d} to {0, 1}. The probability space and measure are defined similarly as 
before. Each u gives rise to a Laplace operator with the following boundary conditions: if the 
—j-th component of the random variable associated to the vertex x has the value one, then 
the edge f decoupled from x by a Dirichlet boundary condition. If the j-th com- 

ponent of the same random variable has value one then the edge [x, x + ej] is decoupled from 
x by a Dirichlet boundary condition. Conversely, those components of the random variable 
which are zero correspond to Kirchhoff boundary conditions. 
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